Conductance in strongly correlated ID systems: Real-Time Dynamics in DMRG 
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A new method to perform linear and finite bias conductance calculations in one dimensional 
systems based on the calculation of real time evolution within the Density Matrix Renormalization 
Group (DMRG) is presented. We consider a system of spinless fermions consisting of an extended 
interacting nanostructure attached to non-interacting leads. Results for the linear and finite bias 
conductance through a seven site structure with weak and strong nearest-neighbor interactions 
are presented. Comparison with exact diagonalization results in the non-interacting limit serve as 
verification of the accuracy of our approach. Our results show that interaction effects lead to an 
energy dependent self energy in the differential conductance. 
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During the past decade improved experimental tech- 
niques have made production of and measurements on 
one-dimensional systems possible X], and hence led to 
an increasing theoretical interest in these systems. How- 
ever, the description of non-equilibrium transport prop- 
erties, like the finite bias conductance of an interacting 
nanostructure attached to leads, is a challenging task. 
For non-interacting particles, the conductance can be 
extracted from the transmission of the single particle 
levels 0, 0, Since the screening of electrons is re- 
duced by reducing the size of structures under investiga- 
tion, electron-electron correlations can no longer be ne- 
glected. Recently several methods to calculate the zero 
bias conductance of strongly interacting nanostructures 
have been developed. One class of approaches consists 
in extracting the conductance from an easier to calcu- 
late equilibrium quantity, e.g. the conductance can be 
extracted from a persistent current calculation Q, from 
phase shifts in NRG calculations Q, or from approxi- 
mative schemes based on the tunneling density of states 
0. Alternatively one can evaluate the Kubo formula 
within Monte-Carlo simulations ||, or from DMRG cal- 
culations 0. In contrast, there are no general methods 
available to get rigorous results for the finite bias conduc- 
tance. While the problem has been formally solved by 
Meir and Wingreen using Keldysh Greens functions [l0( , 
the evaluation of these formulas for interacting systems 
is generally based on approximative schemes. 

In this work we propose a new concept of calculating 
finite bias conductance of nanostructures based on real 
time simulations within the framework of the DMRG 0, 
E1E1E1E1E1E1E3- It provides a unified description 
of strong and weak interactions and works in the linear 
and finite bias regime, as long as finite size effects are 
treated properly. 

In a first approach of real time dynamics within 
DMRG, Cazalilla and Marston integrated the time- 
dependent Schrodinger equation in the Hilbert space ob- 
tained in a finite lattice ground state DMRG calcula- 
tion [lflj . Since this approach does not include the den- 



sity matrix for the time evolved states, its applicability 
is very limited. Luo, Xiang and Wang [1J| improved the 
method by extending the density matrix with the contri- 
butions of the wave function at intermediate time steps. 
Schmitteckert 01 showed that the calculations can be 
considerably improved by replacing the integration of the 
time dependent Schrodinger equation with the evaluation 
of the time evolution operator using a Krylov subspace 
method for matrix exponentials and by using the full fi- 
nite lattice algorithm. 

An alternative approach is based on the wave function 
prediction ^| . There one first calculates an initial state 
with a static DMRG. One iteratively evolves this state 
by combining the wave function prediction with a time 
evolution scheme. In contrast to the above mentioned 
full t-DMRG, one keeps only the wave functions at two 
time steps in each DMRG step. In current implementa- 
tions the time evolution is calculated by approximative 
schemes, like the Trotter decomposition or the 

Runge-Kutta method [lij . In our work, we combined the 
idea of the adaptive DMRG method with direct evalua- 
tion of the time evolution operator via a matrix expo- 
nential as described in Ref. 17] . Therefore our method 
involves no Trotter approximations, the time evolution is 
unitary by construction, and it can be applied to models 
beyond nearest-neighbor hopping. 

The Hamiltonian for the nanostructure attached to 
leads, H = H g + + He is given by 



H s = -fs(c]c i _ 1 +H.c.) + ^figUj 

.7=71+1 j=n 

m— l 
+ V l n ; 

.7=71+1 

h l = Yl -t(c]c j _ 1 +H.c), 

Kj<n,m<j<M 

H c = ~tc (<4 c „-i +4 l c m _ 1 +H.c). 
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Individual sites are labeled according to Fig. Ms = 
to — n is the size of the interacting nanostructure, [ig 
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FIG. 1: Nanostructure attached to leads and schematic den- 
sity profile of the initial wavepacket at T = 0. 



denotes a local external potential, which can be applied 
to the nanostructure, and V is a nearest-neighbor in- 
teraction term inside the nanostructure. The hopping 
elements in the leads, the structure, and coupling of 
the structure to the leads are t, ts, and tc respectively. 
Similar to the approach in 0] we add an external source- 
drain potential Hi = Hsd/2 (j2jZi ™j - Z)ji m n i) > 
to the unperturbed Hamiltonian H and take the ground 
state |^(T = 0)) of H+H\, obtained by a standard finite 
lattice DMRG calculation, as initial state at time T = 
[l7j . In addition we target for the ground state of H. In 
actual calculations the switched external potential was 
smeared out over three lattice sites. 

We then perform a time evolution as described above 
by applyin g th e time evolution operator U = e lHT on 
|*(T = 0)) |23 S3, which leads to flow of the extended 
wave packet through the whole system until it is reflected 
at the hard wall boundaries as described in [l7|. 

The expectation value of the current at each bond and 
every time step is given by 

J jt ,•_!(«) = -^Re{i(^{t)\t. elegit))}. (2) 

Following Refs. [l3ll20| we define the current through the 
nanostructure as an average over the current in the left 
and right contacts to the nanostructure 

J(T) = [J n>n -i(T) + J m , m -i(T)]/2. (3) 

For the calculation of the DC-conductance through the 
nanostructure the time evolution has to be carried out 
for sufficiently long times until a quasi-stationary state is 
reached and the steady state current J can be calculated. 
If the stationary state corresponds to a well-defined ap- 
plied external potential hsd, the differential conductance 
is given by g(fisu) — e 9J(/xsd)/9msd- In the limit of a 
small applied potential, /j,sd — > 0, the linear conductance 
is given by g(fi g ) = eJ(^ g )/^ S D- 

We first consider the transport through a single im- 
purity. The current rises from zero and settles into 
an oscillating quasi-stationary state (Fig. |2Jl ■ After the 
wavepackets have travelled to the boundaries of the sys- 
tem and back to the nanostructure, the current falls back 
to zero and changes sign. The amplitudes of the oscil- 
lations depend on /isD and fi g , and are proportional to 
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FIG. 2: Current through a single impurity with tc = 0.5t at 
nominal half filling N/M = 0.5 obtained from exact numerical 
diagonalization. (a) For different bias voltages (fi g = —0.2). 
(b) For different gate voltages (/^sd = 0.4). System size is 
M = 96. (c,d) For different system sizes: (c) /i 9 = —0.2, and 
(d) fig = 0. (See text for details.) 



the inverse of the system size 1/M. The period of oscil- 
lation strongly depends on the applied potential [Fig. [21 
(a)] but is independent of the gate potential and system 
size [Fig. |2(b), (c)], and is given by T osc = 27rft/|/x SD |). 
This periodic contribution to the current is reminiscent 
of the Josephson contribution in the tunneling Hamilto- 
nian, obtained by gauge transforming the voltage into a 
time dependent coupling tc(T) = tc e lAISr,T . It is present 
even for zero gate potential, but the currents in the left 
and right leads oscillate with opposite phase and cancel in 
the current average Eq. © . After the wavepackets have 
finished one round trip, the current oscillations reappear 
because of the additional phase shift due to the differ- 
ent lengths of the left and right leads [Fig. [21 (b)]. The 
stationary current is given by a straightforward average, 
because the oscillation period T osc is known. In general, 
the density in the leads, and therefore also the current, 
depends on the system size and a finite size analysis has 
to be carried out in order to extract quantitative results 
[Fig. [21(c), see also discussion of Fig. UJ. Only in special 
cases (symmetry, half filled leads, and zero gate poten- 
tial) is the stationary current independent of the system 
size [Pig. 0(d)]. 

Our result for the conductance through a single im- 
purity in Fig. is in excellent quantitative agreement 
with exact diagonalization results already for moderate 
system sizes and DMRG cutoffs. Accurate calculations 
for extended systems with interactions are more diffi- 
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FIG. 3: Current and differential conductance as function of 
applied potential through a single impurity with /i s = and 
half filled leads: N/M = 0.5. Circles (squares) show results 
for tc = 0.5t (0.35t). System size was M = 48 (M = 96) and 
ncut — 200 (400) states were kept in the DMRG. Lines are 
exact diagonalization results for M = 512. 



cult, mainly for two reasons: 1.) The numerical effort 
required for our approach depends crucially on the time 
to reach a quasi-stationary state. For the single impu- 
rity, the quasi-stationary state is reached on a timescale 
proportional to the inverse of the width of the conduc- 
tance resonance, 4th/ Vq, in agreement with the result in 
Ref. [i(J- In general, extended structures with interac- 
tions will take longer to reach a quasi-stationary state, 
and the time evolution has to be carried out to corre- 
spondingly longer times. 2.) In the adaptive t-DMRG, 
the truncation error grows exponentially due to the con- 
tinued application of the wave function projection, and 
causes the sudden onset of an exponentially growing er- 
ror in the calculated time evolution after some time. This 
'runaway' time is strongly dependent on the DMRG cut- 
off, and was first observed in an adaptive t-DMRG study 
of spin transport by Gobert et al.[21|. We observe the 
sudden onset of an exponentially growing error in our 
calculations as well, Fig. but in addition to the de- 
pendence on ncutj the 'runaway' time now also depends 
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FIG. 4: Exponentially growing error in the adaptive t-DMRG 
for large bias voltages and times: Current through a non- 
interacting 7 site nanostructure with tc = 0.5i, is = 0.8t, 
and fig = 0. System size is M = 144 and N/M = 0.5. The 
number of states kept in the DMRG were ncut = 600 (circles) 
and 1000 (squares). Lines are exact diagonalization results. 
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FIG. 5: Differential conductance as a function of bias volt- 
age through a 7 site nanostructure with nearest neighbor 
interaction. Parameters are tc = 0.5t, ts = 0.8t, and 
N/M=0.5. Squares (circles) denote weak (strong) interaction 
with V/ts = 1 (3). Lines are fits to a Lorentzian with an 
energy dependent self energy E = ir]o + ir/ifi 2 . Dashed lines: 
r?i = 0. System size is M = 144 (M = 192) and 600 (800) 
states were kept in the DMRG. 



strongly on /^sd- To avoid these problems one has to re- 
sort to the full t-DMRG [17|, which does not suffer from 
the runaway error. A detailed analysis of the numerics 
of our approach will be published elsewhere [24| . 

In Fig. [5] we show results for the first differential con- 
ductance peak of an interacting 7-site nanostructure. 
Careful analysis of the data shows, that in order to repro- 
duce the line shape accurately, one has to introduce an 
energy dependent self energy. Since the effect is small, we 
approximate it by a correction quadratic in the bias volt- 
age difference fi = /igD — Mpcak- It is important to note 
that for the strongly interacting nanostructure, V/ts = 3, 
the conductance peaks are very well separated. There- 
fore the line shape is not overlapped by the neighboring 
peaks, and the fit is very robust. Performing the same 
analysis for a non-interacting nanostructure with a com- 
parable resonance width, we obtain negligible corrections 
to rji in the self energy, indicating that the change of the 
line shape is due to correlation effects. 

The linear conductance as a function of applied gate 
potential can be calculated in the same manner, if a suffi- 
ciently small applied external potential is used. We study 
the same non-interacting 7-site nanostructure as before 
and use a bias voltage of hsb — 2 • 10~ 4 . For half filled 
leads, the result for the linear conductance calculated 
with a fixed number of fermions, N/M = 0.5, is qualita- 
tively correct, but the conductance peaks are shifted to 
higher energies relative to the expected peak positions at 
the energy levels of the non- interacting system (Fig.|HJ). 
Varying the gate potential \x g increases the charge on the 
nanostructure by unity whenever an energy level of the 
nanostructure moves through the Fermi level [Fig. H3 (b)] . 
The density in the leads varies accordingly [Fig. [fj] (c)]. 
Since the number of fermions in the system is restricted 
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FIG. 6: Transport through a non-interacting 7-site nanos- 
tructure with tc = 0.5t and ts = 0.8t. The energy levels of 
the nanostructure are indicated by dashed vertical lines, (a) 
Linear conductance for different N. The result after applying 
finite size corrections is shown as straight line (see text for 
details), (b) Number of fermions on the 7-site nanostructure. 
(c) Density p = (N - N 0ot )/(M - M s ) in the leads. System 
size is M = 96 and the number of states kept in the DMRG 
is ncut = 400. 



to integer values, direct calculation of the linear conduc- 
tance at constant p is not possible and one must resort 
to interpolation. Using linear interpolation in p(N,/i g ) 
for N = 44 ... 48 yields our final result for the linear con- 
ductance at half filling [Fig. [|J] (a)]. The agreement in 
the peak positions is well within the expected accuracy 
for a 96 site calculation. Our results for the conduc- 
tance through an interacting extended nanostructure are 
presented in Fig.0 The calculation for the weakly inter- 
acting system requires roughly the same numerical effort 
as the non-interacting system. In the strongly interact- 
ing case, where the nanostructure is now in the charge 
density wave regime, the time to reach a quasi-stationary 
state is longer, and a correspondingly larger system size 
was used in the calculation. In both cases we obtain peak 
heights for the central and first conductance resonance to 




FIG. 7: Linear conductance through an interacting 7 site 
system with tc = 0.5t and ts = O.St for weak (squares) and 
strong (circles) interaction. System size is M = 96 (M = 192) 
and 400 (600) states were kept in the DMRG. Finite size 
corrections have been included. Lines are guides to the eye. 



within 1% of the conductance for a single channel. 

We have introduced a new concept of extracting the fi- 
nite bias and linear conductance from real time evolution 
calculations. Very accurate quantitative results are pos- 
sible, as long as finite size effects are taken into account. 
Our results for the linear conductance compare favor- 
ably both in accuracy and computational effort with the 
DMRG evaluation of the Kubo formula |9j . Calculations 
of strongly interacting systems show correlation induced 
corrections to the resonance line shape. 
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